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Abstract 

Expressions  for  the  chemical  potential  of  an  Einstein  solid,  and  of  ideal  Fermi 
and  Bose  gases  in  an  external  one-dimensional  oscillatory  trap,  are  calculated 
by  two  different  methods  and  are  all  found  to  share  the  same  functional  form. 
These  derivations  are  easier  than  traditional  textbook  calculations  for  an  ideal 
gas  in  an  infinite  three-dimensional  square  well.  Furthermore,  the  results 
indicate  some  important  features  of  chemical  potential  that  could  promote 
student  learning  in  an  introductory  course  in  statistical  mechanics  at  the 
undergraduate  level. 

(Some  figures  in  this  article  are  in  colour  only  in  the  electronic  version) 


1.  Introduction 

The  thermodynamic  identity  [1]  describes  differential  changes  in  the  internal  energy  U  of  a 
system  as 

dU  =  TdS  +  (-PdV  +  fidN  +  ■■■),  (1) 

where  T,  S,  P,  V,  /x  and  N  are  the  temperature,  entropy,  pressure,  volume,  chemical  potential 
and  number  of  particles  (assumed  to  be  of  a  single  species  as  opposed  to  a  mixture), 
respectively.  This  expression  assumes  that  the  system  is  close  to  equilibrium,  as  otherwise  it 
will  not  have  a  single  well-defined  temperature,  pressure  and  so  on.  If  the  process  is  reversible, 
one  can  identify  TdS  as  the  heat  added  to  the  system  and  the  sum  of  the  terms  in  parentheses  in 
equation  (1)  as  the  work  done  on  it1 .  If  a  process  is  irreversible,  one  cannot  always  divide  the 
energy  change  into  a  heat  term  and  a  work  term;  an  example  is  kinetic  friction  in  which  the 
mechanical  and  thermal  transfers  are  commingled  at  mesoscopic  length  scales  [2], 

1  Specifically  —  P  dV  is  the  expansion/compression  work  and  //d,V  is  the  diffusive  or  chemical  work.  The  ellipsis 
in  equation  (1)  indicates  that  many  other  kinds  of  work  can  also  occur,  such  as  magnetizing  a  paramagnet  or  charging 
up  a  capacitor. 
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Of  the  various  terms  in  equation  (1),  probably  /i  is  the  least  familiar  to  undergraduate 
students  being  introduced  to  thermodynamics  for  the  first  time  [3,  4].  Part  of  the  reason  may 
be  that  there  are  few  examples  of  analytic  calculations  of  /i  for  a  system  of  particles  other 
than  for  a  classical  ideal  gas  in  a  three-dimensional  box  with  rigid  walls  (typically  computed 
by  first  deducing  the  Sackur-Tetrode  equation  by  approximating  the  translational  partition 
function  as  a  Gaussian  integral,  not  a  trivial  derivation  for  novices). 

Given  that  /i  is  usually  first  introduced  in  thermodynamics  to  describe  diffusion  [5],  i.e. 
variations  in  the  number  of  particles,  a  student  might  reasonably  ask  why  /r  varies  with  T  for 
fixed  N7  In  reply,  one  can  discuss  the  analogous  issue  of  why  P  varies  with  T  for  fixed  V  for  a 
gas.  Part  of  the  difficulty  may  be  one  of  language:  textbooks  state  that  pressure  controls  any 
changes  in  volume  (and  likewise  chemical  potential  controls  any  changes  in  the  number  of 
particles)  which  can  lead  to  the  incorrect  conclusion  that  if  the  number  of  particles  does  not 
change,  then  the  chemical  potential  can  be  ignored. 

In  order  to  improve  student  understanding,  in  this  paper  two  oscillator  models  (a  solid  and 
a  fermionic  or  bosonic  quantum  gas)  are  analysed  for  which  //  is  readily  calculable.  Different 
methods  are  used  to  determine  //  for  these  two  models  and  the  results  are  then  compared  to 
demonstrate  the  underlying  unity  of  the  ideas. 


2.  Chemical  potential  of  an  Einstein  solid 


Einstein’s  model  of  a  solid  [  1  ]  consists  of  a  set  of  N  independent  one-dimensional  oscillators 
(corresponding  to  N/3  atoms  in  three  dimensions)  among  which  a  total  of  q  —  U / hv  integer 
units  of  energy  are  distributed.  (Here  h  is  Planck’s  constant  and  v  is  the  fundamental  frequency 
of  oscillation,  assumed  to  be  the  same  for  all  N  oscillators.)  That  is,  if  one  writes  the  energy  of 
the  ;th  oscillator  (for  i  —  1 , . . . ,  N)  as  e,  =  n,hv  (after  subtracting  off  the  thermodynamically 
uninteresting  zero-point  energy),  where  n,  is  the  quantum  number  of  that  oscillator,  then 
i v 


(2) 


i=i 


In  the  usual  ‘sticks  and  dots’  representations  of  the  possible  microstates  of  this  system,  the 
units  of  energy  are  pictured  as  dots  partitioned  by  (V  —  1  sticks  into  N  groups  representing  the 
oscillators.  The  multiplicity  Q  is  then  equal  to  the  number  of  ways  of  choosing  to  distribute 
the  q  dots  among  the  q  +  N  —  \  total  objects,  which  is  given  by  a  binomial  coefficient 
(q  +  N  -  1)! 


Q  = 


(3) 


q'.(N  -  1)! 

The  Boltzmann  entropy  is  S  —  k  In  Q,  where  k  is  Boltzmann’s  constant.  Assuming  that 
both  q  and  N  are  large,  the  factors  of  unity  in  equation  (3)  can  be  dropped  and  Stirling’s 
approximation,  ln(« !)  ~  n  In/?  —  n  for  sufficiently  large  integer  n,  can  be  invoked  to  obtain 


k  In 


q  +  N\q  (q  +  N 
N 


(4) 


Next,  by  rearranging  equation  (1),  the  standard  statistical  mechanics  definitions  of  temperature 
and  chemical  potential  are  deduced. 


and 


(5a) 

(5b) 
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Figure  1.  Plots  of  the  dimensionless  chemical  potential  as  a  function  of  the  dimensionless 
temperature.  The  upper  curve  is  for  the  ideal  Fermi  gas,  while  the  lower  curve  is  for  either 
the  Einstein  solid  or  the  ideal  Bose  gas. 


where  the  subscripts  indicate  variables  that  are  held  fixed  during  the  partial  differentiation. 
Volume  is  not  an  independent  variable  in  the  Einstein  model,  so  the  V  subscripts  can 
be  dropped  in  equations  (5a)  and  (5b).  Substituting  U  =  qhv  and  equation  (4)  into 
equation  (5a)  leads  to 


N  ehv/kT  _  l’ 

which,  according  to  equation  (2),  equals  the  average  quantum  number  («)  per  oscillator  as  a 
function  of  temperature.  As  T  — »■  0,  this  expression  correctly  predicts  that  (. n }  — »■  0  because 
every  oscillator  is  near  its  ground  state2.  At  the  opposite  limit,  as  T gets  large  (compared  to  the 
characteristic  Einstein  temperature  TE  =  hv/k ),  equation  (6)  predicts  that  U  =  qhv  &  NkT, 
consistent  with  the  equipartition  theorem  (remembering  that  each  one-dimensional  oscillator 
has  one  kinetic  and  one  potential  degree  of  freedom  at  high  temperatures).  Finally,  substituting 
equation  (4)  into  (5b),  performing  the  derivative  and  then  substituting  equation  (6)  into  the 
result  to  eliminate  q  in  favour  of  7’ lead  to  a  compact  expression  for  the  temperature  dependence 
of  the  chemical  potential  for  an  Einstein  solid, 

fi  =  kT  ln(l  —  e_rfi/r).  (7) 

Defining  the  dimensionless  chemical  potential  /i  =  n/ hv  and  temperature  T'  =  T/TE,  then 
equation  (7)  becomes 

fi'=T  ln(l  -  e~1/r)  (8) 

as  plotted  in  figure  1 . 

3.  Chemical  potential  of  a  quantum  gas  in  an  external  parabolic  trap 

Consider  a  set  of  N  noninteracting  particles  confined  in  one  dimension  by  an  applied  harmonic 
potential  so  that  they  translationally  oscillate  at  fundamental  frequency  v.  Their  energy  levels 

2  Technically  T  cannot  be  set  to  zero  because  Stirling’s  approximation  for  q !  was  used  in  the  derivation.  In  practice, 
that  restriction  can  be  ignored  because  if  N  is  of  the  order  of  say  Avogadro’s  number,  then  there  is  a  large  range  of 
values  of  q  for  which  q  / N  1  while  simultaneously  maintaining  q  1 . 
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are  nondegenerate  (ignoring  spin  degeneracy)  and  uniformly  spaced  by  hv,  so  that  the  density 
of  states  g  is  constant, 

g=  7^.  (9) 

hv 

where,  as  in  section  2,  the  particle  energies  e  are  measured  relative  to  the  ground  state, 
eground  state  =  0.  The  total  number  of  oscillators  must  be  equal  to 

n  OO  i  nOO 

N  —  I  £(£)/(£)  de=  —  /  /  (e)  de,  (10) 

Jo  hv  Jo 

where /is  the  occupancy  or  distribution  function.  In  writing  this  as  an  integral  rather  than  as 
a  sum  over  states,  it  is  assumed  that  the  spacing  between  levels  (compared  to  typical  thermal 
energies)  is  small  enough  that  they  form  a  quasi-continuum,  i.e.  hv  kT . 

First  apply  this  formula  to  a  set  of  fermions.  The  Fermi-Dirac  distribution  function  is 


/(£)  = 


1 

e(£— M)/ kT  +  1  ’ 


(ID 


where  /x  is  the  chemical  potential.  Substituting  this  function  into  equation  (10)  and  changing 
variables  to  x  =  (e  —  iT)/kT  give 


kT 

N  =  — 
hv 


l 


°°  dx  kT  „ 


(12) 


As  x  —*■  oo  note  that  ln(l  +  ex)  — »•  ln(eA  )  =  x,  and  hence  the  integral  evaluated  at  the  upper 
limit  is  equal  to  zero.  Therefore  equation  (12)  becomes 


N  =  —  +  —  ln(l  +  e-^kT).  (13) 

hv  hv 

Now  evaluate  this  at  T  —  0,  for  which  /x  =  sF,  the  Fermi  energy.  The  Fermi  energy  is  equal  to 
the  energy  of  the  highest  occupied  level  at  absolute  zero  and  is  therefore  necessarily  positive 
for  our  choice  of  the  zero  of  energy,  so  that  equation  (13)  becomes 

N  —  —  =>•  Sp  =  N  hv  (14) 

hv 

as  expected  because  the  levels  are  nondegenerate.  Substituting  this  result  into  equation  (13) 
leads  to 


sF  =  fi  +  kT  \n{\  +  er^kT), 


(15) 


which  rearranges  into 

/x  =  kT  ln(e7>/7'  -  1), 


(16) 


where  the  Fermi  temperature  is  Tp  =  Sp/k.  Note  that  equation  (16)  correctly  reduces  to 
jjL  —  Sp  in  the  limit3  as  T— >•  0.  On  the  other  hand,  /x  =  0  when  T  =  Tp/ In  2,  in  contrast  to  an 
ideal  Fermi  gas  in  a  three-dimensional  box  with  rigid  walls  for  which  /x  =  0  when  T  ~  Tp 
(cf  figure  7.16  in  Schroeder  [1],  obtained  by  a  nontrivial  numerical  calculation). 

The  dimensionless  chemical  potential  and  temperature  can  be  defined  as  /x'  =  /x / — 
H/Nhv  and  T’  =  T /Tp  —  kT /Nhv,  respectively,  so  that  equation  (16)  becomes 

/x'  =  T  ln(e1/r  -  1)  (17) 

as  is  plotted  in  figure  1.  Defining  \jle  to  be  the  right-hand  side  of  equation  (8)  and  /x'F  to  be 
the  right-hand  side  of  equation  (17),  it  can  be  shown  that  (i!F  —  fi'E  +  1.  In  other  words,  a 


3  Technically  T  again  cannot  actually  equal  zero  because  the  integral  approximation  in  equation  (10)  requires  that 
T  hv/ k.  Nevertheless  if  is  large  enough,  this  requirement  can  be  maintained  even  when  T  <$C  Tp  =  Nhv/k. 
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dimensionless  graph  of  chemical  potential  for  an  ideal  gas  of  fermions  in  one  overall  harmonic 
oscillator  potential  has  the  same  shape  as  that  for  a  collection  of  individual  Einstein  oscillators 
each  in  their  own  potentials.  The  upward  shift  in  the  vertical  intercept  for  the  fermion  case 
is  required  by  the  Pauli  exclusion  principle.  If  one  is  only  interested  in  changes  in  the  value 
of  the  chemical  potential  with  temperature,  as  in  many  chemistry  problems,  then  this  vertical 
shift  is  irrelevant. 

Next,  consider  a  set  of  N  bosons  in  the  same  external  one-dimensional  simple  harmonic 
potential  so  that  they  again  oscillate  at  fundamental  frequency  v.  The  Bose-Einstein 
distribution  function  is 


1 

■f{£)  ~  e(E-M)/A-r  _  j  ’ 


(18) 


where  //  <  eground  state  =  0  is  a  necessary  requirement.  That  is,  the  chemical  potential 
must  be  negative  at  all  nonzero  temperatures,  unlike  for  fermions,  as  otherwise  /  would 
diverge  to  infinity  at  some  (positive)  energy.  However,  it  must  be  true  that  \i—±  0  as  T->  0 
because  otherwise  the  occupancy  of  all  levels  would  be  zero  at  absolute  zero;  in  contrast, 
H  m  £ groUnd  state  permits  all  /V bosons  to  occupy  the  ground  state  at  sufficiently  low  temperatures. 
(Note,  however,  that  Bose-Einstein  condensation  cannot  occur  in  a  one-dimensional  gas  [6].) 
Substituting  equation  (18)  into  equation  (10)  and  again  changing  variables  to  x  =  (e  —  [T)/kT 
gives 

kT  r°°  dx  kT  r 

n=TT  1— T  =  +  ln(e  -  1)]“/*r  (19) 

hv  J _^/ kT  eA  -  1  hv  w 

and  the  integral  evaluated  at  the  upper  limit  is  again  zero4.  Equation  (19)  therefore  becomes 

N  —  —  — —  —  \n(e'fl/kT  -  1),  (20) 

h  v  h  v 

which  rearranges  to 

fi  =  kT  ln(l  -e~Nhv/kT).  (21) 


Introducing  /x'  =  jjt/Nhv  and  T'  =  kT /Nhv  as  for  the  fermion  case,  equation  (21)  becomes 


jT  =  V  ln(l-e~1/r), 


(22) 


which  is  identical  to  equation  (8).  In  fact,  the  similarity  extends  beyond  the  normalized  version 
of  this  expression  to  the  fully  dimensional  forms.  Equations  (7)  and  (21)  are  the  same  if  the 
Einstein  frequency  i >E  is  related  to  the  frequency  of  oscillation  of  the  harmonically  trapped 
bosons  i >b  according  to  vE  =  Ni>b .  Now  consider  typical  experimental  values.  The  Einstein 
temperature  of  copper  [7]  is  240  K  so  that  hvE  =  21  meV.  On  the  other  hand,  the  average 
energy  of  a  particle  in  an  oscillatory  trap  of  effective  spring  constant  k  is  (s)  =  \kL 2  where 
L  is  the  mean  amplitude  of  the  motion  and  thus  represents  the  one-dimensional  volume  of 
the  system  of  noninteracting  particles.  In  particular,  for  the  harmonic  axis  of  a  magnetic  trap 
designed  by  Hess  [8]  for  N  —  1012  hydrogen  atoms  (having  a  mass  of  m  —  1.67  x  10~27  kg), 
(s)  =  2  x  10-27  J  and  L  =  5  cm  so  that 

N  1 2(e) 

NhvB  =  —TiJ  —  =  20meV  (23) 

L  V  m 

which  is  comparable  to  21  meV.  Thus  the  chemical  potential  of  an  Einstein  solid  is  equivalent 
to  that  of  an  ideal  gas  of  bosons  in  an  external  harmonic  potential  both  in  theory  and  in 
practice3 . 


4  Note  that,  since  the  chemical  potential  is  negative,  the  range  of  integration  does  not  encompass  x  =  0  where  the 
integral  diverges. 

5  Notice  from  equation  (23)  that  Nv  is  proportional  to  the  number  density  N/L  of  particles  in  the  trap,  so  that  /;  in 
equations  (16)  and  (21 )  properly  scales  intensively  with  the  temperature  and  density.  In  accordance  with  equipartition, 
(e)  —  kT  at  temperatures  high  enough  that  the  integral  approximation  in  equation  ( 10)  holds. 
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4.  Conclusions 

For  either  an  Einstein  solid  or  a  quantum  gas  in  a  one-dimensional  oscillatory  trap,  the 
dimensionless  chemical  potential  /i'  has  the  same  functional  dependence  on  temperature, 
decreasing  monotonically  with  an  asymptotic  —T'  In  T  dependence  for  large  values  of  the 
dimensionless  temperature  T .  The  calculations  of  /i  presented  in  this  paper  are  more 
straightforward  than  those  for  a  quantum  (or  even  a  classical)  gas  in  a  three-dimensional 
rigid  box  and  thus  could  be  instructively  presented  in  an  introductory  statistical  mechanics 
course  as  a  prelude.  Notably,  the  vertical  intercepts  and  general  shapes  of  the  graphs  of  /x(7) 
are  the  same  for  gases  in  either  a  one-dimensional  harmonic  potential  or  a  three-dimensional 
infinite  square  well  and  have  been  neatly  explained  by  Cook  and  Dickerson  [3].  Furthermore, 
the  fact  that  the  Einstein  solid  and  ideal  Bose  gas  give  the  same  result  when  computed 
by  different  methods  (by  differentiating  the  entropy  in  section  2  and  by  normalizing  the 
distribution  function  in  section  3)  harmonizes  different  concepts  about  chemical  potential  and 
different  models  of  oscillators.  Finally,  presentation  of  the  one-dimensional  Bose  gas  in  a 
simple  harmonic  potential  can  serve  as  a  springboard  for  treatment  of  the  three-dimensional 
case  that  corresponds  to  actual  Bose-Einstein  condensation  experiments  for  dilute  atomic 
gases  [9,  10]. 
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